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Abstract 

The cosmology of the string effective action, including one loop 
string threshold corrections, is analyzed for static compactifications. 
The stability of the minima of a general super symmetry breaking po- 
tential is studied in the presence of radiation. In particular, it is shown 
that the radiation bath makes the minima with negative cosmological 
constant unstable. 
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The gravitational sector of closed string theories [Q] is a model indepen- 
dent prediction that can be used to check the effective theory at low energies. 
In particular, the dilaton has been recently recognized to play an essential 
role in the duality properties of string theory in the presence of time depen- 
dent backgrounds [2-5]. It also plays a distinghished role in the cosmologi- 
cal solutions of the gravitational equations of motion |J , corresponding to a 
scalar-tensor theory of gravity , and providing a time dependence to gauge 
couplings |§. However, at the scale of supersymmetry breaking, a potential 
for the dilaton can develop. In that case, the constancy of gauge couplings 
is related to the stability of the minimum of the scalar potential. 

In a recent paper, Tseytlin and Vafa || analyzed the cosmology of the 
tree level string effective action including the dilaton field. They proposed a 
mechanism, based on the annihilation of winding modes [fT0|| , that allowed the 



radiation era to follow from the stringy phase, and studied the cosmological 
equations of motion in the absence of supersymmetry breaking. They showed 
that the radiation era was an attractor for certain range of initial conditions. 
In this case the dilaton field, and therefore gauge couplings, remain constant 
and general relativity is recovered. 

In this paper we study the cosmology of the string effective action in- 
cluding one loop string threshold corrections. We find that the cosmological 
equations of motion have essentially the same form as those of ref. ||. We 
then introduce a general supersymmetry breaking potential for the dilaton 
and analyze the stability of its minima in the presence of radiation. In partic- 
ular we show that minima with negative cosmological constant are unstable. 

We consider a critical (D = N c + N+1 = 10) heterotic string compactified 
to N = 3 space dimensions in a time dependent compactification background, 
Gmn =diag(g M; ,, e a ^5 mn ). The string effective action in the Einstein frame 



can be written in supergravity form [|Tl|] as 

1 - 



S = f d A x^g y~R - Q\ d,X t d»X> - e g (g, Q^Q 3 - 3) 



+ S m , (1) 



where X; L = (S, T) are scalar fields in chiral supermultiplets, Q = K + ha \ w\ 2 
where K is the Kahler potential and w the superpotential, and S m includes 
the contribution from the rest of (stringy) matter. The S and T fields are 
related to the dilaton and moduli through |T2"| 

ReS = e 3 *- 2 ^ . , 

ReT = e a . ^ > 
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The Kahler potential |T2| , |T~3|] , including one loop string threshold correc- 
tions [14], can be written as [15 



(S, T) = - In (S + S - 3c ln(T + T)) - 3 ln(T + T) 



(3) 



where c = ^ 2f=i ^gs * s the Green-Schwarz anomaly |L3], [L^, [L6|]. In terms 
of the new variable 



Y = S + S - 3c\n(T + T) 



(4) 



the effective action can be written as 



S 



1 



1 (dY s 



4 V Y 



4 



Y 



+ Sn 



(5) 

where we have fixed, for simplicity, Im S'and Im Tto constant values. (For a 
cosmological analysis including a non-trivial Im S, see ref . ]T7[] . ) The action 
(^D can be written in the string frame through a dual invariant |l5j conformal 
redefinition of the four dimensional metric g^ u = Yg^ u [] 



S 



J d 4 xy/^g~Y 



2 2 v y 



4 



y 



(da) 2 -YV(Y,a) 



-S„ 



(6) 

We will study the cosmology associated with this theory, with metric g^ v = 
diag (—1, a 2 (t)5ij), i = 1, ...,iV = 3 , in a similar way as in ||, by defining 
the variables 

Xi(t) = In di(t) 

(7) 



and writing the effective action as 



\nY -J2\nai(t) 



5 = / rftv^; e -* <r £ a, 2 + 2 - 1 (i - f) * 2> J - 2yy(y, a) 



+ <5 m . 

(8) 

1 Note that the action (||) has the structure of a Jordan-Brans-Dicke theory with 
= -1 fl. 
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From now on we will consider static compactification backgrounds (a = 
constant) and study the cosmological evolution of the loop corrected string 
effective action (^j) in the presence of an arbitrary scalar potential and a ther- 

1 dE(X) 

mal bath of stringy modes with total energy E(X) and pressure -P(A) = — — — — — 

JXI OA 

We take the same scale factor in all (N) spatial directions (Aj = A). The 
action can then be written as 



S = f dtyf^Z {e"* [g°° (-NX 2 + y? 2 ) - W(Y) 



where F is the free energy || and W(y) 
motion derived from the action (M) are 



F}, (9) 
2YV(Y). The equations of 



2 - NX 2 = e^E + WiY) 
X - ipX = \ e^P + \YW'(Y) 
(p-NX 2 = \ e^E - \YW\Y) 

along with the energy conservation equation 

E + NXP = . 



(10) 



(11) 



We consider, for large A, an energy E = e x and pressure P = j^E .f\ We 
can write the equations of motion (fLOD as 



x 



■ N 2 
y = — x 
y 2 



-y 

2 y 



-W 2 (z) 



(12) 



z = —z(y + Nx) 



where we have defined new variables y = (p, x = A (a; is the Hubble param- 
eter), z = Y and the potentials W\ and W 2 are given by 



W\(z) = zW'(z) - ±W(z) 
W 2 (z) = zW'(z) + W(z) 



(13) 



2 In the absence of a potential, see eq.([L0[), the solution (p = —NX is an attractor of the 
equations of motion This solution exactly corresponds to a radiation type universe, 
where Y = — Y(ip + NX) = 0, thus leading to constant gravitational and gauge couplings. 
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Since we are interested in general stability conditions we have to analyze 



the non-linear equations of motion (12). The stability of the critical points 



x l of first order differential equations x l = f l {x 3 ) is usually determined from 

df, 

the analysis of the linearized system, x l = M^x 3 , where Ada = — — : L» j 

ox 3 

and in particular from the eigenvalues of the matrix A4ij Q [0|. However, in 



our case the linearized system is not a good approximation, see eq.(|12j), and 
we need to develop new techniques. 

Let us first consider the analysis of critical points and stability of the 
equations ([12]) in the absence of potential (for arbitrary z). They reduce to 

1 2 1 2 

x = --x +xy + —y 

n\ i, < 14 > 

The only critical point of eq.flHD is the origin (x, y) = (0, 0). The eigenvalues 
of the matrix A4 at the critical point are both A = . Therefore the origin 
is a star point but its stability is undetermined fTSI . 



An exact analytical method for studying this type of non-linear homo- 
geneous differential equations is based on the change of variables to polar 

y 

coordinates x = rcosO and y = rsin# so that A = — = tan#. The sta- 
re 

bility can then be studied in terms of the parameter A. The equations of 
motion (0) can be written as 



2 „ NA(A 2 + iV + 2) + A 2 - N 

C0S9 2N(1 + A 2 ) 

(N + A)(N — A 2 ) 



(15) 



6 = r cos 9 



2N(1 + A 2 ) 



It is clear from the second equation that there are fixed directions (9 = 0) 
through the origin for some particular values A = A Q . The critical point is 
attractive (r < 0) or repulsive (r > 0) along these directions depending on 



3 The eigenvalues of this matrix are the slopes of trajectories in phase space and therefore 
determine the stability or instability of those critical points. 
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the sign of y 



A Q = — iV ==>- r = ry 



_ 2N (16) 

A = iviV ==>- r = ry 

Therefore, for y < the origin is stable while for y > it is unstable 
along the fixed directions A Q . However, we are interested in the stability 
of the radiation era, which corresponds to the asymptota y = —Nx. The 
eigenvalue analysis |l9j only gives information about the critical points but 



says nothing about the stability or instability of the fixed directions (A = A Q ). 
We have developed an analytical method for studying the stability of these 
asymptotas, which can be generalized to any non-linear homogeneous system 
of first order differential equations. 

Let us perturb the angle (8 = o + e) to an infinitesimally close direction 
A(e) = tan(0 o + e) ~ A Q + e(l + AA. The angular variation can easily be 
computed as 

■, . . . AN + A{e)){N - A{e) 2 ) N - 2NA Q - ?>A 2 2 , 

g(£)=rSlgn(3;) 2iY(l + W 2 6X ^ +0{€ ) • 

(17) 

An asymptota (A = A Q ) will be stable if an infinitesimally close direction 
in phase space will approach it, i.e. if 9(e) is negative (positive) for positive 
(negative) e. In other words, stable for 6 '(e = 0) < and unstable for 
6 '(0) > (where the prime denotes derivative with respect to e). In our case 
we have 

N — 1 

A =-N ^9 '(0) = -x — — 

/ 1 x (18) 



A = ±VN =}► 6'{0) = -y 1± 



N , 



and therefore the radiation era (A Q = —N) is stable for x > 0, in the region 
between the unstable asymptotas A Q = ±\/N, as shown in the flow diagram 
of Fig.l, i.e. the projection of four dimensional phase space on the plane 
(x, y) []. The differential equations ( p!4|) determine that the trajectories in 
phase space will not escape that region. However, as we can see from our 

4 In ref . J^] the stability of the radiation era was studied only for initial conditions in the 
region x > 0. 
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analysis, initial conditions in other regions will not converge to the radiation 
era. In particular, for x < the radiation asymptota is unstable. This result 
is crucial for understanding the instability of the minima of a scalar potential 
in the presence of radiation. 

We now analyze the cosmological evolution described by the complete 
equations ([12]) . The critical points are f\ 




'W(O) 



W(0) unstable (stable) improper nodes 
unstable saddle points 



and 



y = —Nx 



4 



{ (N - l)z W'{z ) = 2W(z ) 
2 . W(z ) 



{ Vo 



N(N-1) 

-NXn 



(19) 



(20) 



The critical points ( P0|) are the most interesting ones since they correspond 
to the radiation era. For N = 3, the first condition reduces to 



z W'(z ) - W(z ) = 2z 2 V'(z ) = 



(21) 



which gives the extrema of the scalar potential in the Einstein frame, as one 
would naturally expect from dynamical arguments. However, the existence 
of a critical point at the radiation asymptota crucially depends on the sign 
of W(z ). Only for positive or zero cosmological constant (A = V(z )) we 
have a critical point. The solutions of the (cubic) .My— eigenvalue equation, 
for N = 3, 



'W(z ) 



6 



\ 2 + (2W(z ) + z 2 W"(z ))\+4 ] 



W(z ) 



6 



z 2 W"( Zn ) = 



(22) 



determine the stability of the critical points (f20D. A stable center or spiral 
point corresponds to complex eigenvalues, while an unstable saddle point 



5 Note that for W(0) < there are no critical points. This is a general feature of these 
equations. 
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corresponds to real eigenvalues of different sign |nj. Therefore, the con- 
dition of stability of these critical points (which are maxima or minima of 
the scalar potential V(z)) is the existence of complex solutions to the cubic 
equation ([22]) , which gives the condition on the coefficients, and therefore on 
the potential, for N = 3, 

8z 2 W"(z ) - 3W(z ) > . (23) 

This condition can be expressed in terms of the scalar potential V(z) 
as 8z 2 V"(z ) > 3V(z Q ) and interpreted as the well known stability condition 
of extrema of potentials in a curved background, Am 2 > 3 A |2(J, where 
2z 2 V"(z n ) is the mass of the scalar field in the Einstein frame. This 



-T-,2 



ni 



result, obtained from the stability analysis of the gravitational equations of 
motion at the minimum of the potential in the Jordan frame, is in agreement 
with the result obtained in ref. ||20|| from the positivity of the total energy 



functional of the gravitational system. 

So far we have analyzed the stability of the critical points (PD|). Let us 
now analyze the stability of the radiation asymptota in the presence of an 
arbitrary scalar potential. Suppose that we are driven to the minimum of 
the potential W(z ). The differential equations are homogeneous and 
a similar analysis as that of eq. (|TT| ) can be performed. It gives in this case, 
for A a = -N, 

'W-^ + ^w) ■ 

The radiation asymptota is stable (unstable) for large positive (negative) x, 
as in the absence of potential fll%p . The equations of motion ( P"2"|) along the 
radiation asymptota reduce for N = 3 to 

x = -2x 2 + l -W(z ) . (25) 

We have plotted this function in Fig. 2 for different values of W(z Q ). It is 
a typical example of a structurally unstable differential equation [|TjJ. For 
W(z ) = there is just one critical point, and the universe will approach 
it after an infinite time. For W(z ) > there are two critical points, at 

xf = ±J w<y ^ . The positive root is attractive and the other is repulsive. 
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Assuming initial conditions x(0) > (as in the zero cosmo logical constant 



case) the universe will approach the attractor at x 



which corresponds 



to an inflationary universe. For W(z ) < there is no critical point (as we 
know from the previous analysis) and the Hubble parameter decreases until 
it changes sign and the universe starts contracting. 

For x(0) > 0, eq. (|25|) can be integrated to give the Hubble parameter and 
scale factor 



x 



(+) 



,(+) 



^coth2,fe 
6 V 6 



1/2 



(26) 



a;H(t) 




(27) 



(28) 




that correspond to W Q = W(z ) positive, zero and negative respectively. We 
have plotted these functions in Fig. 3. For W Q < 0, the universe will start 
contracting at a time At after the beginning of the radiation era given by 



(29) 



(For a cosmological constant A ~ 10~ 120 in Planck units, we obtain the age 
of the universe, as expected pi]). After the universe has started contracting 
in the presence of a radiation bath, the stability condition (pl|) determines 
that the asymptota A = —N is unstable. From the third equation in (|T2]) 
we deduce that the value of the scalar field Y increases (decreases) for 
y + Nx < (> 0), and escapes from the minimum of the potential. It is 
the radiation bath in an anti-de Sitter background which acts as a source of 
energy for the scalar field to pass the potential barrier. 
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We now apply the previous analysis to a class of supersymmetry breaking 
potentials in string theories, namely those provided by gaugino condensation 
P^] . The existence of a local minimum in the supersymmetry breaking po- 
tential requires at least two gaugino condensates, and the scalar potential 
can be computed from the superpotential 

W (S) = d ie - aiS + d 2 e~ a2S , (30) 

where d{ and an depend on the gauge group of the hidden sector. We will 
not worry about the details but take a generic case which may be good 
phenomenologically [p3 . The scalar potential will have a local minimum 



only for non-trivial imaginary part of the field S. In this case the scalar 
potential can be written as 

YV(Y) = d 2 ie -^ Y (( ai Y + l) 2 - 3) + d 2 2 e- a * Y ((a 2 Y + l) 2 - 3) 

-2rf lC / 2 e- (ai+Q2)y/2 (aia 2 y 2 + (e*i + a 2 )Y - 2) . ' ' ' ' 

where Y is the string loop corrected variable (Q). These potentials have a 
negative cosmological constant at its local minimum, which can be used to 
fix the dilaton and therefore the gauge couplings. We have plotted V(Y) and 
VK(V) in Fig. 4 for particular values of d{ and aj. 

We have solved the non-linear equations of motion ( |1"2"D with Mathematica© 



24j| and shown the phase space projection on the (x, y) plane in Fig. 5 . We 
assumed that the initial conditions are that of radiation at the minimum of 
the potential (y = — 3x, Y = Y Q ). This is a reasonable asumption since the 
scalar potential will begin to be important at a (condensation) scale much 
below the Planck scale and therefore the universe had enough time to con- 
verge to the radiation asymptota. On the other hand, the minimum of the 
potential is stable during the expansion of the universe and therefore the 
scalar field had time to settle at its minimum. As predicted from the general 
analysis (0), during the expansion phase the radiation is stable, but once 
the universe starts contracting (x < 0) it becomes unstable and the solutions 
diverge from the radiation asymptota, thus moving the scalar field away from 
the minimum of the potential, as shown in Fig. 6. Depending on the initial 
conditions (above or below the asymptota) the universe will expand or con- 
tract in finite time, while the scalar field moves towards Y = or Y = oo 
respectively. This is shown in Fig. 7. For a realistic supersymmetry breaking 
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potential based on gaugino condensation, the minimum has a negative cos- 
mological constant of order A ~ m \/2 ~ in Plank units and therefore 

the universe would have collapsed by At ~ 1CT 27 s, much before primordial 
nucleosynthesis. 

In conclusion, we have studied the cosmology of the string effective action 
including one loop string threshold corrections and analyzed the stability of 
the minima of a supersymmetry breaking scalar potential in the presence of 
radiation. For positive or zero cosmological constant the minimum is stable 
and the universe enters an inflationary phase or approaches the radiation 
era, respectively. For negative cosmological constant the minimum becomes 
unstable when the universe starts contracting and the scalar field moves away 
from the minimum of the potential. 
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Figure Captions 

Fig.l Phase space diagram (solid lines) of equations (0), in the projected 
plane (x,y). The stability of the asymptotas, A = ±\/3 (dashed lines) 
and A = —3 (dashed-dotted line), in different regions follows the anal- 
ysis of eq.(|l8|). 

Fig. 2 Phase space (x, x) of eq.(j2"5|) for W > (solid curve), W = (dashed 
curve) and W < (dashed-dotted curve). For W Q < the Hubble 
parameter x will never stop decreasing. 

Fig. 3 The solutions x(t) and a(t) (Hubble parameter and scale factor) of 
eqs.(|26H28l) for W Q > (solid curves), W = (dashed curves) and 
Wq < (dashed-dotted curves). For W Q < the scale factor starts 
contracting after a time interval At given by eq . (|29|) . 

Fig. 4 The scalar potentials V(Y) (solid curve) and W(Y) (dashed curve), 
in the Einstein and Jordan frame respectively, for the hidden gauge 
group SU(4)ixSU(5)2 as determined by c = 0, a\ — 27r 2 , a<i = 16n 2 /5 
and d\ = l/8n 2 e, d,2 = 5/327r 2 e p3| . The local minimum corresponds 
to a negative cosmological constant. Y is given in units of a± while 
W(Y) and V(Y) are in units of d\. 

Fig. 5 Phase space of equations ([12] ) in the presence of the scalar poten- 
tial W(z) of Fig. 4 in the projected plane (x,y), for initial conditions 
very close to the radiation asymptota. As predicted by our analysis, 
the radiation era is stable in an expanding universe but unstable in 
the contracting phase (x < 0) which can only be reached for negative 
cosmological constant. 

Fig. 6 Phase space projection on the (z, x) plane, for initial conditions x = 
1 ± 0.1, y — —3 and z = z Q (very close to the radiation asymptota) in 
the presence of a negative cosmological constant. The minimum of the 
potential is unstable since the value of the scalar field decreases to zero 
(solid curve) or increases to infinity (dashed curve) while the Hubble 
parameter goes to ±oo for y(Q) + 3x(0) > or < respectively. 

Fig. 7 Evolution of the scale factor a(t) corresponding to the initial con- 
ditions of Fig. 6. The universe will expand (solid curve) to infinity or 
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collapse (dashed curve) in finite time after entering the contracting 
phase. 
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